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Abstract. We calculate and resum the nuclear enhanced power corrections 
from the final state partonic scattering in nuclear matter to open charm pro- 
duction and correlations. In p+A reactions, we find that single and double 
inclusive D mesons can be suppressed as much as the neutral pions from the 
dynamical high twist shadowing. Effects of initial state energy loss in p+A 
collisions are also investigated and may lead to significantly weaker transverse 
momentum dependence of the nuclear suppression. 
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1. Introduction 

A good probe of the dense nuclear medium created in relativistic heavy ion collisions 
should be sensitive to its dynamical scales |21 Ejj while it can be cleanly measured 
experimentally and reliably calculated theoretically. Because of color confinement, 
only hard probes can be reliably calculated in the perturbation theory of Quantum 
Chromodynamics. On the other hand, typical dynamical scales of the dense medium 
produced in p+A and A+A reactions are of the order of hundreds of MeV, which 
is much smaller than the scale of a hard probe and therefore non-perturbative. 
Therefore, an ideal probe should be not only "hard" but also sensitive to this 
"soft" physics. Open charm production has the potential to be an ideal probe 
because of two distinctive scales of the open charm meson: the charm quark mass 
{rric ~ 1.5 GeV, a relatively hard scale) and the binding energy (~ Md — rric ~ 
hundreds MeV) [ |2I ■ This contribution summarizes part of our study of nuclear 
effects on D meson production and D-meson-triggered correlations in p+A reactions 
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at RHIC. 



2. Coherent power correction and initial state energy loss 



We calculate the open charm (Z)" and meson) production cross sections to lowest 
order in perturbative QCD, using standard CTEQ6 parton distribution functions 
and fragmentation functions derived by Braaten et al. in heavy quark effective field 
theory [ ^ . The contribution from scattering of the radiatively generated charm is 
included in our calculation El ■ We find that for inclusive single charm production 
c + g ^ c + g and c + c + q{q) dominate over gluon fusion g + g ^ c + c. The 

g + g ^ c + c process controls charm-anticharm pair production. More specifically, 
the cross sections are calculated as follows [j2j: 



da 



NN 



dyid^PT^ 



abed ^ 
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Xa^h 



^|M,,^,,|2 ,(1) 
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The nucleus-induced modification to the cross sections measured in A-\- B col- 
lisions relative to the binary scaled p + p result is identified through the ratio 



R 



(n) 

AB 



da 



AB 



/dyi- ■ ■ dyndpTi ■ ■ ■ dpT„ 



coll\ 
AB/ 



da 



NN 



7^2/1 • ■ -dyndpTi ■ ■ -dpT^ 



(3) 



{-^ab) ™ ^1- ® often evaluated from an optical Glauber model. Here hi can be 
any combination of final state hadrons, for example hi — + /12 = £>^ + . 
The effects of coherence can be cleanly studied in deep inelastic scattering [ 
If the longitudinal momentum transfer becomes small, the probe will interact 
simultaneously with more than one nucleon. In a frame where the longitudinal 
size of the exchange virtual photon is given by 1^ = 1/xP the Lorentz contracted 
nucleon has longitudinal size 27-0/7 — 2r'o/(-P/r7ijv). The critical value for the onset 
of coherence then reads xjy = l/(2romAr) ^0.1 [ 3 . In our work [Ej, we combined 
the computation of the dynamical mass, generated by multiple final state scattering, 
for heavy and light quarks [I^ with the observable effects of high twist modification 
of hadron production in p-l-A reactions Specifically, isolating the small Xb 

dependence in Eqs. ^ and Q into Fab^cd{xb) = 4'b/N{xb)\Mab^cd\^ / Xb, we resum 
all of the high twist nuclear-enhanced corrections in the regime \t\ <C |s|, We 
find 



Fab^cd{xb) 



Fab — >cd I Xb 



(4) 
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Fig. 1. Left panels: coherent power correction effects on inclusive D meson pro- 
duction at two different rapidities, y = 1.25 and y = 2.5, in d+Au collisions at 
^/sNN = 200 GeV. Right panels: the effect of power corrections and energy loss on 
D meson production. Bottom panels show the centrality dependence of the nuclear 
effects for y = 2.5 only. 

In Eq. I0J Cd = 1 (9/4) for quarks (gluons), respectively. When ^ we recover 
the result of [0. (Note that can be the mass of the charm quark.) 

The left panel of Figure 1 shows the nuclear modification RpA^^^ (pti) fo^' 
the single inclusive Z?*^ + production for two different rapidities yi = 1.25 and 
yi = 2.5. The nuclear effect is similar to that for light pions The centrality 
dependence of high twist shadowing is also shown. Up to rapidity yi = 2.5, the 
effect of power corrections is relatively small. Preliminary PHENIX data on 
the nuclear modification for muons, presumably coming from heavy quark decays, 
is shown for comparison. The nuclear effect appears to be larger than the one 
predicted by power corrections. Similar results for i?^^ ^(ptdPTs) in 

{D^ + £)+) — (D° + D~) correlations as a function of the rapidity gap, y2 — yi, and 
centrality are given in the left panels of Figure 2. 

Additional nuclear attenuation may come from the energy loss of incoming par- 
tons in cold nuclear matter. Initial state medium induced radiation, when followed 
by a hard scattering in the nuclear environment, must be clarified in the pQCD 
formalism prior to precise quantitative applications in cold nuclear matter. Thus, 
we are motivated to carry out only a phenomenological study of AE effects by the 
known incoherent Bertsch-Gunion type radiation, which scales with the path length 
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through the nucleus, L [[7] 

^-^ ^.J . AEi^^^C.^E^m-^, (5) 
dLud''k±_ TT"^ kj_(k - qi)j_ tt Xg ^^qcd 

where fi — L/Xia the mean number of scatterings. If the incoming partons lose 
a fraction of their energy, the effective rapidity of the collision is shifted in the 
direction of the propagation of the nucleus (backwards for RHIC kinematics). For 
minimum bias collisions, we use a rapidity shift of Ay — 0.25. This is not incompat- 
ible with the nuclear modification of light hadrons [ 6 and the rapidity asymmetry 
observed in low energy p+A reactions (^/s = 20 GeV), where effects of coherence are 
negligible [E]- The variation of Ay with centrality is given by the nuclear thickness 
Tb(6). The right panels of Figures 1 and 2 show the consequence of such an effective 
implementation of the energy loss. Firstly, the nuclear suppression has a weak pT 
dependence in both single inclusive open charm production and charm anticharm 

correlations. Secondly, the effect on R^pA ^(ptuPTs) 

can be large, compared to power corrections for transverse momenta larger than 
few GeV. In particular, the right panel of Figure 1 shows suppression much more 
compatible with the preliminary PHENIX measurements. A constant rapidity shift 
may fail if we go to either high px or forward rapidity since we do expect energy 
dependence of the fractional energy loss [E]- This points to the necessity of careful 
studies of energy loss effects in cold nuclear matter. 

We note, finally, that it is a characteristic of both power corrections and initial 
state energy loss that they reduce the rate of the hard scattering processes, thus 
similarly affecting single and double inclusive hadron production. 



3. Conclusions 



In summary, we have presented the first calculation of coherent resummed QCD 
power corrections to open charm production and open charm triggered back-to 
back correlations [E]. We find that the nuclear suppression is qualitatively similar 
to that for light hadrons Initial state energy loss effects have also been phe- 
nomenologically incorporated and shown to yield a different transverse momentum 
behavior of the nuclear suppression. We expect our theoretical results will guide the 
experimental determination of the relevant many body QCD effects that control the 
dynamics of forward rapidity hadron production in p-l-A reactions [HI 1^. They will 
also provide the much needed baseline for precision heavy quark QGP tomography [ 
^ at forward rapidities. 
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Fig. 2. Nuclear attenuation from power corrections (left panels) and power correc- 
tions and energy loss (right panels) in correlated back-to-back D and anti D meson 
production in d-|-Au at RHIC. The trigger D meson rapidities are yi — 1.25 and 
yi = 2.5 and the associated anticharm mesons are measured at 1/2 = 0. Centrality 
dependence is shown for yi = 2.5 only. 
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